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I.  INTRODUCTION 


An  experimental  program  on  the  three-dimensional  response  of  concrete 
cubes  (Ref.  1)  was  recently  completed  at  New  Mexico  State  University  (NMSU). 
The  results  were  presented  as  plots  of  stress  versus  strain  for  each  concrete 
specimen  tested.  No  in-depth  evaluation  of  the  quality  and  consistency  of  the 
data  was  performed  by  NMSU.  The  New  Mexico  Engineering  Research  Institute 
(NMERI)  subsequently  evaluated  the  data  obtained  by  NMSU.  This  report  pre¬ 
sents  revised  results  based  on  that  evaluation. 

In  the  NMSU  experiments,  each  concrete  cube  specimen  was  subjected  to  a 
prescribed  stress  path  test  until  the  specimen  failed.  Each  type  of  stress 
path  test  was  performed  on  two  or  three  concrete  specimens.  Because  multiple 
runs  had  been  performed  for  each  path,  it  was  possible  for  NMERI  to  inspect 
the  data  for  apparent  discrepancies  and  baseline  shifts.  The  data  were  then 
adjusted,  and  composite  stress-strain  curves  were  developed.  These  curves 
should  prove  to  be  invaluable  for  the  development  and  evaluation  of  constitu¬ 
tive  models  for  plain  concrete. 

In  Sections  II  and  III  of  this  report,  the  stress  paths  on  which  the 
NMSU  experimental  program  was  based  are  described,  and  experimental  values  for 
limit  states  are  presented.  The  paths  and  limit  points  are  given  in  the  two 
stress-invariant  spaces  that  are  used  for  constitutive  modeling. 

The  actual  corrections  made  by  NMERI  in  the  stress-strain  data  and  the 
reasoning  on  which  the  changes  were  based  are  summarized  in  Section  IV.  Plots 
of  both  the  original  and  the  adjusted  data  are  given  so  that  a  reader  who 
interprets  the  results  differently  from  the  authors  can  apply  independent 
judgment.  Composite  stress-strain  curves  are  given  for  each  stress  path  in 
Section  V. 


1.  Traina,  L.  A.,  Experimental  Stress-Strain  Behavior  of  a  Low  Strength 
Concrete  Under  Multiaxial  States  of  Stress,  AfWL-TR-82-92,  Air  force 
Weapons  Laboratory,  Kirtland  Air  Force  Base,  New  Mexico,  December  1982. 


Invariants  of  inelastic  strain  are  used  in  most  constitutive  models. 

These  invariants  can  also  be  used  to  provide  an  alternate  check  on  the  ade¬ 
quacy  of  the  composite  stress-strai n  curves.  The  check  consists  of  illus¬ 
trating  each  path  in  a  strain-invariant  space.  If  the  strain-invariant  curves 
obtained  by  using  the  composite  stress-strain  data  provide  a  reasonable  repre¬ 
sentation  of  the  strain-invariant  curves  obtained  by  using  corrected  stress- 
strain  data,  the  indication  is  that  the  composite  stress-strain  data  points 
have  been  chosen  in  a  consistent  manner.  The  results  of  this  check  are  dis¬ 
cussed  in  Section  VI.  Limit  states  in  terms  of  strain  invariants  are  also 
given.  The  conclusions  drawn  as  a  result  of  the  evaluation  program  are  pre¬ 
sented  in  Section  VII. 
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II.  STRESS  PATHS  AMD  LIMIT  STATES  INVOLVING  FIRST  AND  SECOND  INVARIANTS 


In  most  of  the  engineering  literature  concerned  with  the  testing  of 
concrete,  first  and  second  invariants  of  stress  are  used.  These  invariants 
are  the  mean  pressure  (P)  and  the  second  invariant  of  the  stress  deviator, 
(*'"37],  defined  as  follows: 


where  Oj,  a2 ,  and  o3  denote  principal  stresses  positive  in  tension.  These 
parameters  are  used  to  define  the  stress  paths  that  were  employed  in  the  NMSU 
experimental  program. 

The  tests  were  divided  into  three  sets.  Each  of  the  first  set  of  eight 
stress  path  tests  was  performed  on  two  specimens.  The  set  consisted  primarily 
of  conventiona1  uniaxial,  biaxial,  tri axial,  and  shear  tests  with  all  stress 
components  in  compression.  The  second  set  of  nine  tests  involved  more  compli¬ 
cated  paths  in  the  compressive  regime;  consequently,  each  type  of  test  in  this 
set  was  performed  on  three  specimens.  The  last  set  consisted  of  ten  tests  in 
which  one  component  of  stress  was  in  the  tensile  region  for  at  least  some 
portion  of  each  path.  Again,  each  type  of  test  in  this  set  was  performed  on 
three  specimens. 

To  perform  each  experiment,  it  was  necessary  to  give  the  stress  paths 
in  incremental  form.  This  is  also  a  convenient  way  to  describe  the  paths 
analytically.  Tables  1,  2,  and  3  provide  a  complete  description  of  the  paths 
in  terms  of  a  positive  increment  in  stress.  As.  Each  path  is  identified  with 
a  number  and  the  set  to  which  it  belongs,  and  each  segment  of  the  path  is 
described  briefly.  Terms  are  defined  in  the  tables  as  needed.  In  particular, 
f^  and  f£  denote  the  limit  stress  in  uniaxial  compression  and  biaxial  stress. 
The  corresponding  limit  stress  in  uniaxial  tension  is  f£. 
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TABLE  1.  SET  1:  STRESS  PATHS  IN  COMPRESSION 


Pa 

th 

Description 

1 

■  1 

1 

1  Uniaxial  stress  to  failure 

aoj.  =  ao2  =  0;  ao3  =  -as  Limit  stress  f '  =  (-o3) 

c  m  q  x 


2  Uniaxial  stress  to  failure.  This  and  path  1  give  an  initial  test 
for  anisotropy. 

uOj  =  -  As ;  AO  2  =  -  Ao  3  =  0  . 

3  Biaxial  stress  to  failure  with  oJoi  =  1.0;  o3  =  0 

Ao,  =  -As;  ao2  =  -As;  ao3  =  0  Limit  stress  fj  =  (-a,)  =  (-0-) 

1  *  J  b  v  l'max  ^  max 

! 

4  Hydrostatic  loading  to  C;  uniaxial  loading  to  failure 

(i)  uO ;  =  A 02  =  AO3  =  -<aS 

to  the  point  ax  =  o2  =  o3  =  -f^/3 
( i  i )  ^0 ,  -  - -s ;  ^0  )  -  ao  i  -  0 

5  Hydrostatic  loading  to  2/3  f^;  shear  loading  to  failure 

( i )  Ao  1  =  Ao 2  =  Ao3  *  -As 

to  the  point  o1  =  o2  =  03  =  -2/3  f^ 

1 

(ii)  Ao j  =  -As;  A o2  =  As;  ao3  =  0 

6  3-D  stress  to  failure  (subject  to  extensometer  capability)  with 
qJox=  1.0;  o3/ oL  =  0.1 

A01  *  -as;  ao2  =  -as;  uo3  =  -0.1  As 
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TABLE  2.  SET  2:  ADDITIONAL  STRESS  PATHS  IN  COMPRESSION 


Path 

_ 

Description  ! 

1 

(i)  Load 

hydrostati cal ly 

to  03=  o2=  03 

= 

•0.1  f'; 

ao  2  =  = 

AO  3 

=  -uS 

i 

t 

(ii)  Load 

uniaxially  to  a 

*«f;(«i!-o 

.85 

f C ’  ^  =  °3  =  -0.1  f') 

i 

£0  l  = 

-  aS  j  ^2  =  ^03 

=  0 

|  (iii)  Unload  to  01=02=03=  -0.1  P;  &o  1  =  -s;  ao2  =  =  0 

i  (iv)  Load  biaxially  in  xi-x^  plane  to  failure  with 
1  02/ °i  =  1.0;  03  =  -0-1  f ' 

I 

j  <^0  j  =  a0 2  ~  “ as i  A03  =  0 

2  I  (i)  Load  hydrostatically  to  03  =  o2  =  03=  -0.1  f'; 

! 

J  a0 |  =  ao 2  ~  ao 3  -  - uS 

1 

:  (ii)  Load  biaxially  in  x  ;-x2  plane  to  f  ai  lure  witn  =  1.0; 

03  =  -0.1  f ' 

(Check  on  stress-induced  anisotropy) 

AO  |  -  aO ^  ~  “  aS  j  aUg  =  0 

3  |  (i)  Load  hydrostatically  to  01  =  o2  =  03  =  -0.1  fl; 

a 

!  Ao 1  -  a02  =  aCJ  3  =  -as 

(ii)  8i axi al  stress  in  xi-x2  plane  to  0  =  u 
l°i  =  o2  =  -0.85  f£;  03  =  -0.1  j 

ao  ^  -  aO  2  =  -  As  j  ao  3  =  0 

(iii)  Unload  to  ot  =  o2  =  03  =  -0.1  f^;  o;  =  ao2  =  us;  AU3  =  0 
(iv)  Load  uniaxially  in  xi-direction  to  failure 

aO 1  =  -  as j  ao  2  =  *^03  =  0 
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TA8LE  2.  CONTINUED 


Description 


Load  hydrostatically  to  0!  =  o2  =  o3  =  -0.1  f^; 

iio  j  =  Ao2  =  Ao3  =  -  AS 

Load  uni  axially  in  x^-direction  to  failure  (test  of  stress 

induced  anisotropy) 

i^a  i  -  As;  Ao2  =  Aa3  =  0 

Load  hydrostatically  to  u3  =  o2  =  a3  =  -f^ 

Aa1  =  La 2  -  Aa3  =  -As 

Unload  hydrostatically  to  o-  =  o2  =  o3  =  -2/3  f^ 

LxO  ^  ~  AO  2  s  Ao  3  =  AS 

Load  in  shear  to  failure 
lo i  =  -As;  ao2  =  as;  ao3  =  0 

Load  hydrostatically  to  a.  -  a2  *  o.  »  -2  f 


oa  I  =  uo2  =  Aa3  =  -AS 


Unload  hydrostatically  to  =  o2  =  u3  =  -0.1  f^ 

SO|  =  ao2  =  Ao3  =  As 

Load  bi axially  to  failure 
Ao3  =  Ao2  =  -As;  Ao3  =  0 

Load  hydrostatical  ly  to  <jl  -  a2  =  a3  -  -0.1 

Aa  i  =  Aa2  =  Aa3  =  -  As 

Load  uniaxially  to  oI  =  -0.75  f^;  o2  =  u3  =  0 
Acsi  =  -as;  Ao2  =  Ac3  =  0 

Load  to  a  biaxial  state  =  o2  =  -0.75  f^;  o3  =  0 

Ao ^  =  uo3  =  0;  Ao2  =  -as 

Load  biaxially  to  failure 
aoi  =  Aa2  =  -As;  A a3  =  0 


TABLE 


CONCLUDED 


Path 


Description 


(1) 

Load 

hydrostatical ly 

r+ 

O 

Q 

°2  = 

o3  =  -0.1 

Ao  1  = 

=  A02  =  A03  =  -AS 

(11) 

Load 

uni  axial ly  to  oL 

=  -0. 

.35 

f 

C 

Aa1  = 

=  -as;  ao2  =  Ao3 

=  0 

(iii) 

Trace  a  path  to  01  = 

-0.1  1 

:  1  , 

c* 

o2  ■ 

=  -0.35  f ' 

o3  = 

-0.1  f1 

c 

1  ■ 

!  As;  A o2  =  -As; 

Ao  3  = 

0 

(iv) 

Load 

uni  axially  to  fa 

i  lure 

Ao  j  = 

:  ao3  =  0;  Ao2  = 

-AS 

(i) 

Load 

hydrostatical ly 

to  Oi 

= 

o2  = 

03  =  -2  f 

AOj  =  A o2  3  Ao3  =  -As 


(iii) 

(iv) 


(a)  uoi  =  as;  a o2  =  -as;  ao3  =  0  to  ox  =  -f^; 

02  =  -3  f o3  =  -2  f ' 

(b)  aoj  =  -As;  ao2  =  As;  ao3  =  0  to  oj  =  o2  =  o3  =  -2 

Unload  hydrostatically  to  ox  =  o2  =  o3  =  -0.1 
Aa^  =  Ao2  =  ua3  =  As 

Load  bi axi ally  to  failure 
Ao^  =  ao2  =  -As;  Ao3  =  0 


in 


TABLE  3.  SET  3:  STRESS  PATHS  INVOLVING  TENSION 


M 

Description 

— 

1 

Uniaxial  stress  to  failure  o:  =  f£;  o2  =  a3  = 

Ao;  =  As;  Afl2  =  AO3  =  0 

0 

l 

i 

* 

Uniaxial  stress  to  failure  oj  =  o2  =  0;  o3  = 

A<Ji  =  A  o2  =  0;  Ao3  =  As 

(test  for  anisotropy) 

1 

3 

Biaxial  stress  to  failure 

Ao  1  =  As,  Aa2  =  -0.5  As;  Aa3  =  0 

i 

l 

4 

Shear  to  failure  oL  =  f^;  a2  =  -f^;  a3  =  0 

A01  =  As;  Aa2  =  -As;  Aa3  =  0 

5 

! 

1 

, 

(i)  Uniaxial  stress  to  oj  =  0.5  fj.;  o2  -  03 
Ao j  =  As;  Ao2  =  Ao3  =  0 

(ii)  Shear  to  failure 

Ao ^  =  As;  Ao  _  -  -As;  Ao3  =  1 

*  0 

! 

!  6 

i 

1 

(i)  Shear  to  03  =  0.5  f^;  o2  =  -0.5  fj;  03 

Ao ^  =  As;  Ao2  =  -As;  Ao3  =  0 

(ii)  Uniaxial  stress  to  failure 

Ao3  =  As;  Ao2  =  0;  Ao3  =  0 

=  0 

7 

Curved  path  to  failure 

(i)  Ao  2  =  As  Ao  2  =  -0.2  As  A 03  =  0 

(to  01  3  0.5  f£) 

1 

(ii)  Ao 2  =  As  Ao2  =  -0.4  As  Ao3  =  0 

(to  failure) 

8 

Curved  path  to  failure 

(i)  A02  =  As  Ao2  =  -0.3  As  Ao3  =  0 

(to  0i  =  0.5  f^) 

(ii)  A02  3  As  Ao2  =  -0.5  As  Ao3  =  0 

(to  failure) 

TABLE  3.  CONCLUDED 


Description 

(i)  Initial  compressive  state  a2  =  a3 

Aoj  =  0;  Aa2  =  Ao3  =  -As 

(ii)  Uniaxial  stress  to  failure 
Ac^  =  As;  Ao2  =  Ao3  =  0 

(i)  Initial  compressive  state  a3  =  -fj. 

Ao |  =  Aa2  =  0;  Aa3  =  -as 


( i  i )  Shear  to  f ai lure 


The  paths  are  also  shown  in  Figure  1  for  the  /J^-P  plane.  As  these  plots 
demonstrate  rather  vividly,  the  paths  cover  a  large  spectrum  of  possible 
load  so  the  resulting  stress-deformat  ion  relations  should  provide  a  comprehen¬ 
sive  test  for  a  constitutive  model. 

Of  particular  significance  is  the  limit  state;  that  is,  the  state  at 
which  an  incremental  increase  in  deformation  ceases  to  yield  a  corresponding 
incremental  increase  in  stress.  Historically,  this  state  has  been  called 
failure,  but  structural  failure  will  occur  only  if  the  structure  is  statically 
determinant  and  a  limit  state  is  reached.  For  many  materials,  this  limit 
state  is  described  by  a  line  in  the  /J‘-P  plane.  The  limit  points  for  all 

paths  considered  in  this  study  are  shown  in  Figure  2.  This  plot  demonstrates 
that  the  use  of  a  straight  line  for  a  limit  surface  is  a  poor  fit  for  these 
data. 
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(b)  Set  2. 


Figure  1.  Continued. 
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III.  LIMIT  STATES  AND  STRESS  PATHS  INVOLVING  FIRST  AND  THIRD  INVARIANTS 

Lade  (Ref.  2)  suggested  that  a  better  fit  to  limit  data  would  be  obtained 
if  a  third  invariant  of  stress  were  used  instead  of  the  second  invariant. 

This  suggestion  has  been  investigated  (Ref.  3)  for  a  third  invariant,  L, 
defined  by 


L  = 


~  a 


s 


1/3 


where  a$  is  a  parameter  that  can  be  interpreted  as  the  shear  strength  of  the 
material  under  zero  mean  pressure.  Limit  points  in  the  L-P  space  are  shown  in 
Figure  3.  It  is  evident  that  a  straight  line  in  this  space  is  a  much  better 
fit  to  these  data  than  is  a  straight  line  in  the  /jJ-P  space. 


These  results  suggest  that  for  the  formulation  and  evaluation  of  con¬ 
stitutive  models,  the  use  of  the  L-P  space  may  be  preferable  to  the  use  of  the 
/IF-P  space.  Thus,  for  the  sake  of  completeness,  stress  paths  in  the  L-P 
space  are  shown  in  Figure  4  for  each  of  the  three  sets  of  tests.  The  hydro¬ 
static  compression  path  is  included  as  a  reference  curve  because  L  is  not 
i  ’dependent  ;■  ?.  !  po:.is  involving  compression  lie  n  toe  small  sector 

between  the  hydrostatic  compression  curve  and  the  limit  line.  The  advantage 
of  this  presentation  is  that  it  shows  the  similarities  of  many  of  the  conven¬ 
tional  testing  paths  used  for  concrete.  A  disadvantage  is  that  some  resolu¬ 
tion  is  lost  in  comparison  to  that  obtained  in  plots  in  the  space  of  */T^-P  or 

in  plots  in  the  space  of  two  principal  stresses. 


2.  Lade,  P.  V.,  "Three-Parameter  Failure  Criterion  for  Concrete,"  Journal 
of  the  Engineering  Mechanics  Division,  ASE,  108:850-863,  1982. 

3.  Schreyer,  H.  L.,  “A  Third-Invariant  Plasticity  Theory  for  Frictional 
Materials,"  Journal  of  Structural  Mechanics,  June  1983. 
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P,  MPa 


Figure  3.  Limit  states  in  L-P  plane  for  all  paths 
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Figure  4.  Concluded. 


IV.  ADJUSTMENTS  TO  STRESS-STRAI N  DATA 


PRELIMINARY  COMMENTS 

The  data  obtained  from  Reference  I  were  digitized  for  processing  and 
plotting  on  a  Tektronix  mi ni processor  system.  These  data  are  given  in  Appen¬ 
dix  A.  Multiple  runs  for  a  given  stress  path  are  shown  on  one  plot.  The 
stresses  are  positive  in  compression  (negative  in  tension);  the  strains  are 
positive  in  contraction  (negative  in  extension).  The  data  have  been  plotted 
so  that  the  scale  for  contractive  strains  appears  along  the  bottom,  and  curves 
of  contractive  strain  start  at  the  lower  left  corner.  The  scale  for  extensive 
strains  appears  along  the  top,  and  curves  of  extension  start  at  the  lower 
right  corner.  The  use  of  two  scales  yields  a  higher  resolution  of  the  strain 
data  than  is  available  with  a  single  scale. 

The  same  data  are  given  in  plots  of  pressure  versus  volumetric  strain, 
ev,  where 

ev  =  -(ej  +  e2  +  e3] 

I ;  t  o  :h;i  s!gn  c.  ant*  on,  f -e  ■  .  ':etri  ;  ; :•  e  for  como  :;c ‘.  on 

and  negative  for  dilatation.  These  plots  also  appear  in  Appendix  A. 

From  a  survey  of  all  the  uniaxial  stress  path  data,  which  included  four 
concrete  specimens  from  set  1  and  an  additional  four  specimens  from  the  con¬ 
crete  batches  used  in  set  2,  the  elastic  constants  were  determined.  The 
Young's  modulus  suggested  by  these  data  was  21  GPa,  and  the  Poisson's  ratio 
was  about  0.25.  These  values  were  used  when  the  need  for  data  adjustments  on 
each  plot  was  considered.  Adjustments  to  the  data  were  based  only  on  the 
linear  elastic  regime  of  the  data.  Linear  elasticity  and  the  aforement ioned 
elastic  constants  were  used  to  determine  the  strains  from  the  given  stresses. 
These  calculated  strains  were  then  compared  with  strains  obtained  experimen¬ 
tally,  and  a  baseline-type  shift  was  made  in  the  experimental  strains  for  data 
that  were  significantly  different  from  the  calculated  strains.  Some  subjec¬ 
tivity  was  involved  in  these  adjustments,  so  all  changes  were  recorded.  The 
reasoning  behind  the  alterations  is  explained  in  the  succeeding  paragraphs. 
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ADJUSTMENTS  TO  SET  1 


Some  anomalies  and  missing  data  points  were  found  in  several  of  the 
sets  of  strain  data.  The  first  correction  of  the  data  was  made  to  compensate 
for  mechanical  difficulties  attendant  upon  the  measurement  of  small  strains. 
For  this  purpose  the  elastic  portion  of  each  set  of  data  was  compared  to  the 
calculated  theoretical  strains,  and  a  strain  correction  was  added  to  or  sub¬ 
tracted  from  all  of  the  strains  (except  the  initial  zero  strain). 

To  correct  for  anomalous  data  found  for  strain  in  the  1-direction  of 
path  1,  specimen  3-6,  was  set  equal  to  c2  because  this  specimen  had  been 
uniaxially  loaded  in  the  o3  direction,  and  therefore  £ l  should  equal  For 
the  missing  seven  data  points  for  of  path  4,  specimen  D-2,  the  values  of  e3 
were  used  because  this  test  theoretically  produces  strains  of  £2  =  £3*  anti  the 
available  data  points  from  the  other  run  confirmed  these  values.  This  data 
set  and  tne  data  set  for  path  1,  specimen  3-6,  remain  inconsistent.  For  five 
other  specimens,  one  or  two  data  points  were  missing  from  one  axis  of  princi¬ 
pal  stress  at  the  limit  state.  For  these  points,  the  given  stress  data  were 
used  to  find  the  estimated  volumetric  strain,  and  the  volumetric  strain  versus 
lain  pressure  curve.;  we-e  axt-*aoolateJ.  T-.e  :..vc  k '.own  strains  ex¬ 
tracted,  which  left  a  relatively  reliable  value  for  the  missing  strain.  The 
data  corrected  in  this  manner  were  path  3,  specimen  B-3,  e3;  path  3,  speci¬ 
men  C-3,  £3;  path  5,  specimen  D-5,  a3;  path  6,  specimen  C-l,  e3;  and  path  3, 
specimen  D-4,  e3. 

All  the  corrections  are  summarized  in  Table  4. 

ADJUSTMENTS  TO  SET  2  (NONSTANDARD  COMPRESSION  DATA) 

Corrections  to  the  data  for  set  2  consisted  only  of  adjusting  the  ini¬ 
tial  strain  data  to  make  them  consistent  with  linear  elasticity.  The  correc¬ 
tions  are  listed  for  each  path  and  concrete  specimen  in  Table  5.  The  incon¬ 
sistencies  in  the  data  that  made  the  corrections  necessary  were  caused  primar¬ 
ily  by  the  same  mechanical  difficulties  in  measuring  strain  as  were  encoun¬ 
tered  in  the  data  for  set  1. 


TABLE  4.  MODIFICATIONS  TO  SET  1 


Path  Specimen  Adjustment  to  strain  measurement,  ue 

^1  e2  ^3 


1 

C-6 

-70 

1 

B-6 

el  =  e2 

0 

0 

2 

C-2 

0 

0 

Vv'l 

2 

S-2 

0 

0 

3 

C-3 

o 

°  ! 

3 

1  J 

B-3 

0 

0 

0 

I 

!  4 

D-2 

-40 

+300 

0 

A- 2 

o 

+160 

-20 

i 

(n  =  2,3) 

(n  =  2) 

l 

+300 

( n  V,  o\ 

'  v.'  |  | ' 

1 

I 

» 

\  ^  s  / 

V"  - 

-700 
(h  ~  4) 

5 

A-5 

o 

0 

0 

5 

D-5 

0 

+170 

6 

C-l 

-50 

0 

6 

8-5 

-50 

0 

7 

B-4 

0 

HhmPIB 

0 

V"  -  c) 

+250 
(n  >  2) 

7 

C-4 

0 

0 

+150 

8 

D-4 

-80 

-200 

3 

A-5 

-150 

( «  -  ON 

-30 

+50 

m 

■  I X 

Note:  n  represents  the  number  of  the  data  point  in  a  run 
where  n  =  1  corresponds  to  e  =  0  and  o  =  0. 
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Note:  n  represents  the  number  of  the  data  point  in  a  run  where 
n  =  1  corresponds  to  e  =  0  and  a  =  0. 
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CONCLUDED. 


Path 

Specimen 

n 

Adjustment  to  strain  measurement,  ue  j 

£l 

£2 

_ 

£3 

8 

A-3 

=  2 

+120 

+30 

+  150 

=  3 

+  110 

+30  ! 

+  105 

=  4 

+  190 

+  100  ! 

-50 

>  4 

+290 

+30  ; 

-110 

8 

B-3 

=  2 

0 

+50 

=  3 

+  160 

0 

+35 

=  4 

+40 

-30 

>  4 

■ 

+40 

,  | 

-150 

8 

C-3 

!  =  2 

+70 

+60 

+  150 

1  =  3 

+  125 

+  70 

+  135 

=  4 

+80 

+70 

+  170 

>  4 
;  ^ 

'  +70 

i 

+  70 

9 

3-4 

=  2 

+170 

+130 

!  +100 

=  3 

+260 

+210 

+30 

1 

=  4 

+380 

+230 

+  140 

>  4 

+400 

+  130 

o 

C-4 

I  =  2 

+100 

;  +110 

=  3 

—  4, 

j 

r  L  1  ^ 

t  _  j 

j 

>  4 

!  +50 

+200 

9 

i  D-4 

j _ 

L  0 

+200 

! 

f 


Note: 


n  represents  the  number  of  the  data  point  in  a  run  where 
n  =  1  corresponds  to  e  =  0  and  o  =  0. 
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ADJUSTMENTS  TO  SET  3 


Typically,  very  small  strains  are  recorded  for  concrete  in  tension  (on 
the  order  of  100  ye  at  failure,  compared  to  several  thousand  ye  at  compression 
failure).  The  small  strains  are  difficult  to  measure  accurately  with  the 
extensometers  that  were  used  in  these  tests.  The  tensile  data,  therefore,  are 
less  reliable  than  the  compression  data.  For  this  reason,  no  attempt  was  made 
to  shift  anomalous  data.  Instead,  where  data  ootained  for  one  of  the  three 
specimens  were  inconsistent  with  data  obtained  for  the  other  two  specimens, 
the  inconsistent  data  were  neglected.  These  deletions  are  noted  in  the 
strain-stress  plots  in  Appendix  B,  where  all  the  data  are  snown;  the  incon¬ 
sistent  sets  are  marked  with  an  x.  The  deletions  were  based  on  the  strain- 
stress  plots  as  well  as  on  the  plots  of  pressure  versus  volumetric  strain  and 
plots  involving  strain  invariants. 

REVISED  PLOTS 

After  the  adjustments  had  been  made,  the  data  shown  in  Appendix  A  were 
plotted  again.  The  results  are  shown  in  Appendix  8.  Thus,  the  reader  can 

observe  the  effects  of  any  adjustment  by  co-mar’ r a  corresponding  *'i cur'es . 

Most  of  the  tests  for  a  given  path  had  good  repeatabi 1 i ty.  As 
expected,  the  concrete  material  was  not  ideally  isotropic.  This  is  noticeable 
in  the  difference  between  the  uniaxial  tests  in  which  o3  was  the  uniaxial 
direction  (set  1,  path  1)  and  those  in  which  was  the  uniaxial  direction 
(set  1,  path  2).  The  compression  path  for  which  repeatability  was  probably 
least  successful  was  path  9  of  set  2.  This  path  started  with  a  large  hydro¬ 
static  load  followed  by  a  shear  cycle,  hydrostatic  unload,  and  then  biaxial 
loading  to  failure.  The  loading  and  unloading  did  not  produce  completely 
repeatable  strains. 

For  the  P-e7  graphs,  large  dilatation  is  evident  in  the  uniaxial 
(set  1,  paths  1  and  2),  biaxial  (set  1,  path  3),  and  triaxial  (set  1,  path  8) 
test  data  as  well  as  in  most  of  the  data  paths  for  set  2.  The  large 


31 


dilatation  is  particular  to  this  low-strength  concrete;  higher  strength  con¬ 
cretes  exhibit  much  less  dilatation.  (For  example,  see  Reference  4.)  Even 
though  the  triaxial  data  for  set  1,  path  8,  exhibited  a  fair  amount  of  dilata¬ 
tion,  the  triaxial  data  for  set  1,  path  6,  exhibited  almost  none.  Another 
interesting  comparison  can  be  made  of  the  path  4  data  for  sets  1  and  2.  In 
both  cases,  the  specimens  were  subjected  to  a  hydrostatic  load  and  then  a 
uniaxial  load  to  failure.  Set  1,  path  4,  had  a  larger  hydrostatic  load  and 
exhibits  almost  no  dilatation,  while  set  2,  path  4,  with  the  small  hydrostatic 
load,  has  a  large  amount  of  dilatation.  The  behavior  of  the  concrete  speci¬ 
mens  depended  upon  the  amounts  of  hydrostatic  pressure  applied  in  these  two 
stress  path  tests.  The  material  behavior  changed  from  a  uni  axially  loaded 
response  to  something  more  like  a  triaxially  loaded  response. 


4.  Green,  S.  J.,  and  Swanson,  S.  R.,  Static  Constitutive  Relations  for 

Concrete,  AFV1L-TR-72-244,  Air  Force" Weapons  Laboratory,  Kirtlahd  Ai>  Force 
8ase,  New  Mexico,  1973. 


32 


V.  COMPOSITE  STRESS-STRAIN  CURVES 


The  adjusted  strain-stress  data  are  given  in  Appendix  3.  Eacn  stress 
path  test  was  executed  on  two  or  three  different  concrete  specimens,  and  the 
data  for  all  specimens  have  been  plotted  on  one  graph  for  each  stress  path. 
From  these  sets  of  data,  average,  or  composite,  strai n-stress  curves  were 
developed  for  each  stress  path.  The  repeatability  of  the  test  data  was  con¬ 
sidered.  When  the  results  of  one  test  were  very  different  from  those  of  tne 
other  two  tests  in  a  series,  the  anomalous  test  data  were  eliminated  from  the 
development  of  the  composite  curve.  In  the  nonstandard  compression  test  ser¬ 
ies  (set  2),  one  of  the  sets  of  strain-stress  curves  sometimes  differed  only 
slightly  from  the  other  two  sets  of  curves  for  a  particular  stress  path.  In 
such  cases,  the  two  'similar  sets  of  data  were  given  more  consideration,  and 
were  weighted  more  heavily,  than  was  the  slightly  anomalous  curve.  The  more 
inconsistent  data  occurred  in  the  tension  test  series. 


The  composite  curves  were  overlaid  on  the  experimental  data  curves  as 
shown  in  Appendix  S.  The  curves  are  presented  separately  for  stress  versus 
strain  in  Fioures  5,  6,  and  7  and  for  pressure  versus  volumetric  strain  •  ■> 
::,gures  8,  9,  ana  10.  Scalier  pare  are  given  ia  indicate  fr.e  cc^ncc  „r  me 
experimental  data.  These  curves,  which  show  both  the  deformation  and  limit 
features  of  the  material,  are  generally  most  useful  in  the  development  of 
constitutive  models. 
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Set  1  Path  6  (composite) 


Figure  5.  Continued. 


Figure  6.  Continued. 


Figure  6.  Continued 


Strai n 


n  ?^<x 


Set  3  Path  2  (composite) 


Strai  n 

(a)  Paths  1  and  2. 


Figure  7.  Composite  stress-strain  curves,  set  3 
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(b)  Paths  3  and  4. 
Figure  7.  Continued. 
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(c)  Paths  5  and  6. 
Figure  7.  Continued. 
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Figure  7.  Concluded. 
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Path  1  (composite)  Set  1  Path  2  (composite) 


pressure-volumetric  strain  curves,  set  1. 


Set  2  Path  6 
'  (composite) 


Set  3  Path  1  (composite)  Set  3  Path  2  (composite) 


Figure  10.  Composite  press'. re-volunetric  strain  curves,  set  3. 


VI.  DEFORMATION  PATHS  AND  LIMIT  STATES  INVOLVING  STRAIN  PATH  INVARIANTS 

Strain  path  invariants  comprise  an  important  part  of  the  theory  of 
plasticity  if  such  a  model  is  used  to  predict  the  hardening  and  softening 
behavior  of  a  material.  To  adequately  predict  the  limit  state  in  terms  of 
strain,  a  relationship  involving  strain  invariants  is  needed.  No  completely 
reliable  relationship  has  been  defined,  but  approximate  relationships  based  on 
some  form  of  the  strain  invariants  have  been  used  with  adequate  results. 
Because  the  strain  invariants  are  central  to  algorithms  used  in  the  attempt  to 
predict  the  behavior  of  materials,  three  strain  path  invariants  are  considered 
here.  These  invariants  are  defined  as  follows: 

K  °/tr  (d?' ) 

e(  ./  (tr  deiddeidy'2 

e3  =  J  | det  dei | 

1  :  #•  di  ;  4 

where  e'  is  the  inelastic  strain  ,e'  =  e"  -  a",,  e  ~  is  the  inelastic  strain 

f  0 

deviator,  e  is  the  total  (measured)  strain,  and  e  is  the  elastic  strain. 

Plots  of  e^  versus  e^  and  of  e]  versus  e^  for  the  adjusted  stress-strain 
data  are  given  in  Appendix  C.  Corresponding  curves  based  on  the  composite 
stress-strain  data  are  overlaid  on  the  plots  in  the  appendix  and  are  shown 
separately  in  Figures  11  through  16.  The  correlation  between  the  overlaid  or 
composite  curve  and  the  curves  based  on  adjusted  data  represents  a  measure  of 
how  well  the  composite  curves  were  chosen  in  the  first  place.  Overall,  the 
results  shown  in  Appendix  C  attest  to  the  reasonability  of  the  composite 


curves . 


Path  1  (composite)  Set  1  Path  2  (composite) 


■  2  Path  5  (composite) 


(b)  Paths  5  through  8. 
Figure  12.  Continued. 


2  Path  9  (composite) 


5  (composite)  Set  1  Path  6  (composite) 


Concluded 


:  2  Path  6  (composite) 


3  Path  2  (composite) 


3  Path  5  (composite) 


Set  3  Path  10  (composite) 


16.  Concluded. 


The  strain  invariant  plots  of  versus  show  distinct  characteris¬ 
tics  for  the  various  stress  paths.  Most  of  the  compression  paths  indicate  a 
fair  amount  of  negative  inelastic  volumetric  strain.  Two  exceptions  are  found 
in  paths  4  and  6  of  set  1.  Path  4  denotes  a  hydrostatic  load  followed  by  a 
uniaxial  load  to  failure.  Path  6  represents  a  triaxial  stress  path  with 
=  o2  and  a3  =  O.loi.  Both  paths  exhibit  large  inelastic  shear  (or  devia- 
toric)  strains  but  relatively  small  inelastic  volumetric  strains. 

Limit  states  in  terms  of  these  invariants  have  also  been  determined. 

Figure  17  shows  limit  states  in  the  e1  versus  e1  space;  corresponding  points 
i  i  s  v 

are  shown  in  the  e3  versus  ev  space  in  Figure  18.  There  appears  to  be  no 
simple  relation  in  terms  of  the  first  and  second  invariants  to  define  the 
limit  state,  but  the  results  shown  in  Figure  13  in  terms  of  the  first  and 
third  invariants  indicate  a  modest  improvement.  With  the  use  of  all  three 
invariants  it  may  be  possible  to  find  a  simple  strain-based  criterion  that 
defines  the  limit  state.  It  is  particularly  important  that  this  criterion  be 
accurate  for  the  uniaxial  stress  test  because  this  test  is  often  used  as  a 
standard  for  evaluating  material  properties. 
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Figure  17.  Limit 


VII.  CONCLUSIONS 


Except  in  a  few  isolated  cases,  the  quality  of  the  experimental  data  is 
remarkably  good.  No  deviations  from  expected  concrete  behavior  were  observed 
for  any  of  the  paths,  and  strain  responses  in  terms  of  prescribed  stress 
appear  to  be  smoothly  varying  functions.  These  results  should  provide  a  good 
base  for  evaluating  the  application  of  constitutive  models  to  weak  concrete. 

Although  these  data  cover  a  significant  range  of  concrete  behavior,  it 
must  be  emphasized  that  no  data  of  a  comparable  comprehensive  nature  exist  for 
other  regimes  that  involve  strain  softening,  general  cycling  paths,  and  multi¬ 
dimensional  tensile  paths.  The  apparatus  used  in  the  NMSU  tests  could  be  used 
to  provide  the  experimental  capability  for  obtaining  data  in  these  regimes 
with  a  minor  modification  to  provide  displacement-controlled  heads.  Such  data 
are  necessary  if  constitutive  models  are  to  be  verified  for  applications  that 
involve  large  inelastic  deformations. 

A  linear  relationship  between  a  form  of  the  third  invariant  and  the 
mean  normal  pressure  represents  the  limit  state  very  well  in  the  L-P  space.. 
Current  1  j,  limit  states  a.-  e  nor.  represented  wail  o.y  a  s i  —p  1  e  re’atipns.rp  in 

which  strain  invariants  are  used.  However,  there  are  indications  that  alter¬ 
nate  combinations  of  invariants  will  be  successful.  Additional  work  in  this 
area  could  lead  to  the  development  of  plasticity  algorithms  that  are  as  useful 
for  frictional  materials  as  for  metals. 
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APPENDIX  A 


STRESS-STRAIN  CURVES 


This  appendix  contains  the  original  stress-strain  data  for  the  compres¬ 
sion  stress  paths  of  sets  1  and  2.  The  tension  data  are  not  included  because 
no  adjustments  were  made  to  them.  The  original  tension  data,  which  are  the 
same  as  the  adjusted  data,  are  presented  in  Appendix  B.  The  stress-strain 
( a- e)  plots  in  Appendix  A  are  followed  by  mean  normal  pressure-volumetric 
strain  (P-ey )  plots. 
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o-£  Curves,  Set  1  78 
o-e  Curves,  Set  2  82 
P-ey  Curves,  Set  1  87 
P-ey  Curves,  Set  2  89 
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APPENDIX  B 


ADJUSTED  STRESS-STRAIN  CURVES  WITH  COMPOSITE  PLOTS 


Adjusted  stress-strain  curves  with  composite  curves  overlaid  are  pre¬ 
sented  in  this  appendix.  Curves  for  tension  specimens  eliminated  from  consid¬ 
eration  in  the  construction  of  the  composite  curves  are  marked  with  an  X.  The 
mean  normal  pressure  versus  volumetric  strain  curves  also  appear  in  this 
appendix. 
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Adjusted  P-ev  Curves  with  Composite  Curves,  Set  1  108 

Adjusted  P-e^  Curves  with  Composite  Curves,  Set  2  110 

Adjusted  P-ev  Curves  with  Composite  Curves,  Set  3  113 
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LOADING  PATHS  AND  COMPOSITE  CURVES  IN  TERMS  OF 
INELASTIC  STRAIN  INVARIANTS 


This  appendix  contains  all  of  the  paths  plotted  in  two  strain-invariant 
spaces.  The  composite  curves  constructed  from  the  cr-e  paths  and  plotted  in 
the  strain-invariant  spaces  are  superimposed  on  their  respective  experimental 
curves. 
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Adjusted  e^  versus  e^  with  Composite  Curves,  Set  3  123 

Adjusted  e^  versus  e 3  with  Composite  Curves,  Set  1  126 

Adjusted  e^  versus  e3  with  Composite  Curves,  Set  2  123 

Adjusted  e^  versus  e3  with  Composite  Curves,  Set  3  131 
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